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Abstract 

The discrete single-particle spectra in both the Fermi and Dirac sea have been calculated by the 
imaginary time step (ITS) method for the Schrodinger-like equation after avoiding the "tsunami" 
of the Dirac sea, i.e., the diving behavior of the single-particle level into the Dirac sea in the direct 
application of the ITS method for the Dirac equation. It is found that by the transform from the 
Dirac equation to the Schrodinger-like equation, the single-particle spectra, which extend from the 
positive to the negative infinity, can be separately obtained by the ITS evolution in either the Fermi 
sea or the Dirac sea. Identical results with those in the conventional shooting method have been 
obtained via the ITS evolution for the equivalent Schrodinger-like equation, which demonstrates 
the feasibility, practicality and reliability of the present algorithm and dispels the doubts on the 
ITS method in the relativistic system. 

PACS numbers: 24.10.Jv, 21.60.-n, 21.10.Pc, 03.65.Ge 
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The imaginary time step (ITS) method which searches for the direction of the steepest 
descent of energy and follows it in iterative steps until the local minimum on the energy 
surface is reached, is an effective approach to solve the non-relativisticproblems in coordinate 
space and has achieved lots of success in the conventional mean field |2| . As noted already in 
Ref. since only the lowest eigenstate is reached after several exponential transformations, 
the method is applied to the single-particle Hamiltonians of which the spectrum is bounded 
from below. 

From the physics point of view, a lot of novel aspects of nuclear structure and entirely 
unexpected features have been found in the "exotic nuclei" . The extreme neutron 

richness of exotic nuclei and the physics related to the low density in the tails of matter 
distributions have not only attracted attention in nuclear physics and astrophysics but also 
provided a challenge to solve complex many-body problem in coordinate space. 

As one of the best candidates for the description of exotic nuclei, the Relativistic Mean 
Field approach l7(| has achieved lots of success in describing many nuclear phenomena during 



the past years [8l. I9I. llO] 



Subsequently, the Relativistic Continuum Hartree-Bogoliubov (RCHB) theory [Hj pro- 
vides a fully self-consistent treatment of pairing correlations in the presence of the continuum 
and thus a reliable description of nuclei far away from the line of /^-stability, which well re- 
produced the halo phenomena in 11 Li. A natural extension of the RCHB theory is the 
exploration of deformed exotic nuclei and to pin down the existence of deformed halo. Ef- 
forts along this line have been made in the past several years. Due to the difficulty in solving 
the coupled channel differential equations for deformed system in coordinate space [12I ]. the 
method for expansion in Woods-Saxon basis was proposed 131 ]. which becomes very time 
consuming for the heavy system. 

For the ITS method, the extension from spherical to deformed systems is straightforward. 
However, whether the ITS method can be applied for the relativistic system still remains 
an interesting question due to the existence of the Dirac sea. Although, by implementing 
the effective Hamiltonian deduced for the upper component to the existing ITS routine for 
the Skyrme-Hartree-Fock approach, efforts have been made to get the single-particle levels 
in the Fermi sea, e.g., in Ref. [ijj], further works are still necessary to justify its feasibility 
and reliability, and to understand the influence of the Dirac sea explicitly. 

In this paper, the ITS method will be applied to solve the Dirac equation in coordinate 
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space. The disaster of the Dirac sea confronted in the direct application of the ITS method 
for the Dirac equation will be analyzed and discussed in detail. Reliable and practical recipe 
to avoid the disaster will be provided by examining the relation between the Dirac equation 
and its "equivalent" Schrodinger-like equation. 

The main idea and detailed formalism of the ITS method can be found in Ref. The 
ITS evolution of the single-particle wave functions reads 

|^ n+1) ) = (l-A^)|^ ri) >, A = At/h, (1) 

where the exponential series is truncated at the first order as in Ref. j^j, h is the single- 
particle Hamiltonian and {<Pa } is a set of orthogonal single-particle wave functions at time 
t n = nAt. Orthonormalizing {\E f i n+1 ' ) } thus obtained generates a new set of orthogonal wave 
functions {v ? i n+1 ^} which will be used for the next ITS evolution. The total energy of the 
system will decrease as the time evolves, until it finds a local minimum on the energy surface. 
The corresponding wave functions {(p a } will simultaneously converge to the eigenstates of 
the single-particle Hamiltonian h. 

The ITS method has been successfully implemented in the conventional mean field ap- 
proach. For the relativistic case, the main issue is to solve the Dirac equation, 

{a ■ p + (3[M + S(r)} + V(r)} <p a = (e a + M)<p a , (2) 

with the scalar and vector potentials S(r) and V(r). For clarity and simplicity, V(r) ±S(r) 



are assumed to be the spherical Woods-Saxon potentials similar as in Ref. 13[, then the 
Dirac spinor takes the form, 

1 ( l F a (r)y a (r) \ 

<Pa = ~ \ Xi(?a), (3) 

r \G a (r)a r y a (r) j 2 

with the isospinor xi(Qa), an d the spherical spinor y a (r). Throughout this paper, the single- 
particle state in either Fermi or Dirac sea is labeled by the quantum numbers {n a , l a , n a ,j a } 
of the upper component, where K a = ±|j a + l/2| for l a = j a ±l/2. The radial Dirac equation 
for F a and G a can thus be obtained as 



( V + S -# + ^ ^ 

dr r 

,1 + ^ V-S-2M > 
\ dr r / 





(4) 
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FIG. 1: (color online) Disaster for solving the Dirac equation with scalar and vector Woods-Saxon 
potentials in the ITS method. Taking 16 O as an example, the diving behaviors into the Dirac sea 
during the ITS evolution are shown for the neutron single-particle level syi-, which is evolved in a 
spherical box with R = 20 fm, and mesh sizes dr = 0.5 fm (short-dashed line), 0.2 fm (dash-dotted 
line), and 0.1 fm (dashed line) respectively. The shaded areas are the continuum in the Fermi and 
Dirac sea. 

in which the 2x2 matrix operator on the l.h.s. is the single-particle Hamiltonian h to be 
used in Eq. (CQ) for the ITS evolution. 

For the initial orthogonal wave functions {ip^}, the upper components are taken as spher- 
ical Bessel function and the lower components are assumed to be zero. The first (second)- 
order differential operator in the single-particle Hamiltonian is calculated by the seven(nine)- 
points formula as in Ref. Q]. The ITS evolution is carried out in coordinate space within 
a spherical box [0,R] with the mesh size dr = 0.1 ~ 0.5 fm and a typical time step 
At = 10~ 26 s. 

Taking the Woods-Saxon potentials for V(r) ±S(r) corresponding to 16 as an example, 
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the ITS evolutions of the lowest s\/2 orbital for neutron are shown in Fig. HJ It is found that 
although the initial state is set in the Fermi sea, the single-particle energy will inevitably 
dive into the Dirac sea as the time evolves. This is not surprising. Because the ITS method 
searches for the states with the lowest energy in the whole spectra, the single-particle energy 
should dive to the negative infinity corresponding to the wave function with infinite nodes. 
In practice, however, as the node number is limited by the given box and mesh sizes, one 
can always achieve convergence for the ITS evolution. The convergent single-particle energy 
will dive deeper in the Dirac sea for smaller mesh size. 

One may wonder whether the ITS evolution is able to produce any s\/2 orbital outside 
the Dirac sea. Attempts have been made but a dilemma has been found. As indicated in 
Fig. [TJ if one takes a large mesh size, the lowest sy2 orbital may lie shallow in the Dirac 
sea, and the whole spectrum is sparse. Taking a box with R = 10 fm and dr = 0.5 fm as an 
example, a survivor 2s 1/2 is found in the discrete region of the Dirac sea, and is compared 
with the "exact" results obtained by the shooting method 10| with R = 20 fm and dr = 0.1 
fm in Fig. El Its single-particle energy and radial wave functions are in good agreement 
with the "exact" ones, which shows that at least some physical solutions in the Dirac sea 
are accessible by the direct ITS evolution for the Dirac equation. However, it is at a cost 
of precision, and it is gloomy to go beyond the Dirac sea. To calculate more single-particle 
states, one should take smaller mesh size or larger box, which determines the number of 
orthogonal wave functions. However, the whole spectrum of S1/2 orbitals becomes denser 
and dives deeper, thus it is incurable to go beyond the Dirac sea. 

Therefore, this is indeed a great disaster to perform the ITS evolution directly for the 
Dirac equation, since physical solutions are completely submerged by the mighty flood with 
negative energy states, or the "tsunami" of the Dirac sea. 

In order to avoid the "tsunami", one possibility is to perform the ITS evolution for the 
Schrodinger-like equation instead of the Dirac equation. Naively, such possibility might be 
doubted as the Schrodinger-like equation is deduced from the Dirac equation, and they are 
equivalent in the sense that all the solutions of the Dirac equation satisfy the corresponding 
Schrodinger-like equation. However, they are not equivalent in the sense discussed below. 

From Eq. (j4j), the relation between the upper and lower components can be obtained as, 




(5) 
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FIG. 2: (color online) Neutron radial wave functions and single-particle energy for the survivor 
1s\j2 in the disaster of the Dirac sea. The ITS evolution for the Dirac equation is performed in 
a spherical box with R = 10 fm and dr = 0.5 fm, while the "exact" results are obtained by the 
shooting method with R = 20 fm and dr = 0.1 fm. 

with the effective mass M + = M + (S — V + s a )/2, and the corresponding Schrodinger-like 
equation for the upper component is hF a = e a F a , with the effective Hamiltonian, 



* _ Id 2 1 dM + d 
~ ~2M~dr~ 2 + 2M1 dr dr + 



1 dM + K a | 1 K a (K a + l) 



2M 2 dr r 2M, r 2 



(6) 



Now the difference between the Dirac equation and the corresponding Schrodinger-like 

equation is clear. For the Dirac equation, the single-particle spectra extend from the positive 

to the negative infinity. While for the Schrodinger-like equation, although the effective 

Hamiltonian in Eq. ([6]) which contains the single-particle energy is quite sophisticated, it 

P 2 

can be reduced as h = — - + U(e a ). Therefore, with the trial single-particle energy 

2M + {e a ) 

e a in the Fermi sea, the positive definite effective mass will ensure that e a > U(s a ). As 
a result, the single-particle spectrum in the Fermi sea can be separated from that in the 
Dirac sea, and then the ITS evolution could be applied for the Schrodinger-like equation to 
obtain the solutions in the Fermi sea. Furthermore, in order to obtain the solutions in the 
Dirac sea, one could apply the similar technique to perform the ITS evolution for the charge 
conjugated Schrodinger-like equation. Therefore, it is the transform from the Dirac equation 
to the Schrodinger-like equation that allows the ITS evolution for the single-particle spectra 
in either the Fermi sea or the Dirac sea, and thus provides a " life buoy" to avoid the disaster 
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FIG. 3: (color online) Neutron radial wave functions and single-particle energies for lp spin doublets 
in the Fermi sea (upper panels) and 2p spin doublets in the Dirac sea (lower panels) obtained by 
the ITS evolution for the Schrodinger-like equation with R = 20 fm and dr = 0.1 fm. 

of the Dirac sea. 

Using the effective Hamiltonian in Eq. ([6]), the upper component is obtained by the ITS 
evolution for the Schrodinger-like equation. The corresponding lower component is obtained 
by Eq. (JSJ). Then the Dirac spinor including both the upper and lower components are 
orthonormalized for the next iteration. 

As examples, the neutron radial wave functions and single-particle energies for lp spin 
doublets in the Fermi sea and 2p spin doublets in the Dirac sea obtained by the ITS evolution 
for the Schrodinger-like equation with R = 20 fm and dr = 0.1 fm are respectively given in 
the upper and lower panels in Fig. [3j Both the wave functions and single-particle energies 
are identical with the "exact" results obtained by the shooting method. 

Similarly, all the discrete single-particle spectra identical with the "exact" ones for neu- 
tron of 16 O in both the Fermi and Dirac sea are obtained as shown in Fig. HI The spin 
doublets are presented in pairs, with the left (right) one is for k > (< 0). It should be 
noted that the single-particle states are always labeled by the quantum numbers of the upper 
component, different from those in Ref. [ljj]. Therefore, it is clearly shown that the disaster 
of the Dirac sea can be avoided via the ITS evolution for the equivalent Schrodinger-like 
equation. 

In conclusion, the discrete single-particle spectra in both the Fermi and Dirac sea have 
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FIG. 4: (color online) Neutron discrete single-particle spectra of 16 in the Fermi and Dirac sea 
obtained after avoiding the disaster of the Dirac sea via the ITS evolution for the Schrodinger-like 
equation. For each pair of the spin doublets, the left (right) is for k > (< 0). 

been calculated by the ITS evolution for the Schrodinger-like equation after avoiding the 
"tsunami" of the Dirac sea. 

The ITS method, which has been widely used in the conventional mean field approaches, 
will inevitably meet a serious challenge when it is directly applied for the Dirac equation. 
All the single-particle levels will dive into the Dirac sea during the ITS evolution, and the 
physical solutions of the Dirac equation are submerged by the "tsunami" of the Dirac sea. 
Fortunately, the Schrodinger-like equation deduced from the Dirac equation provides the life 
buoy to survive the "tsunami". For a long time, the Schrodinger-like equation was generally 
regarded to be equivalent to the Dirac equation. However, by the transform from the Dirac 
equation to the Schrodinger-like equation, the single-particle spectra, which extend from 
the positive to the negative infinity, can be separately obtained by the ITS evolution in 
either the Fermi sea or the Dirac sea. The ITS evolution for the Schrodinger-like or charge 
conjugated Schrodinger-like equation will provide the solution in the Fermi or Dirac sea, 
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which are identical with the "exact" ones obtained by the shooting method. 

The successful implementation here fully demonstrated the feasibility, practicality and 
reliability for the ITS method in the relativistic system. 
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